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symmetric media equation M. Equation (1) yields the 
two limits 
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The dielectric constant of a conductor-insulator mixture 
shows a pronounced maximum above the critical volume con- 
centration. Further experimental evidence is presented as well 
as a theoretical consideration based on a phenomenological 
equation. Explicit expressions are given for the position of 
the maximum in terms of scaling parameters and the (com- 
plex) conductances of the conductor and insulator. In order to 
fit some of the data, a volume fraction dependent expression 
for the conductivity of the more highly conductive component 
is introduced. 

The ac and dc conductivities of resistor and resistor- 
capacitor (RC) networks and continuum conductor- 
insulator composites have been extensively studied for 
many years. In systems where there is a very sharp 
change (metal-insulator transition or MIT) in the dc 
conductivity at a critical volume fraction or percolation 
threshold denoted by <j> c , the most successful models, for 
both the dc and ac properties, have proved to be ap- 
proaches using percolation theory Early work concen- 
trated on the dc properties, but since it was realized that 
the percolation threshold is a critical point, and that the 
percolation equations could be arrived at from a scaling 
relation, several papers ( Q and the references therein), 
reporting experimental results on the ac conductivity 
have appeared. Review articles, containing the theory 
and some experimental results, on the complex ac con- 
ductivity and other properties of binary metal-insulator 
systems include 0]- Q . 

I n @ ) U an d r| the following equation was introduced 
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which gives a phenomenological relationship between 
ac,cri, and cm- They are, respectively, the conductivi- 
ties of the conducting and insulating component and the 
mixture of the two components. Note that all three quan- 
tities ac, aj, and <tm can be real or complex numbers in 
Eq.(l). The conducting volume fraction <fi ranges be- 
tween and 1 with <j) = characterizing the pure insula- 
tor substance (ctm = <jj) and = 1 the pure conductor 
substance (cm =gc)- The critical volume fraction, or 
percolation threshold is denoted by <p c , where a transi- 
tion from an essentially insulating to an essentially con- 
ducting medium takes place, and A = (1 — <j) c ) /<f> c . For 
s — t — 1 the equation is equivalent to the Bruggeman 
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which characterize the exponents s and t. Note that Eqs. 
(2) and (3) are the normalized percolation equations. At 
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up to higher order terms in aj/ac- 

Thus, in the crossover region where 4> lies between 1 
(<Ji/ac) 1/{s+t) and C + (a//a c ) 1/(s+t) we have [|- 
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When using the above equation to analyze ac systems 
the complex conductivities (a x — a xr — iu;e e rx , where u) 
is the angular frequency, e Q the permittivity of free space 
and e rx the real relative dielectric constant of the com- 
ponent) must be inserted in Eq. (jl|) . Below (j> c the leading 
term yields the imaginary conductivity or real dielectric 
constant, while the next order term gives the real con- 
ductivity. In turn, above <f> c the real conductivity is the 
leading term and the real dielectric constant appears in 
the next order term, upon which the present paper fo- 
cuses. In many instances the approximation ac = a cr 
and 07 = — iuje e r is made @- [@. In this case, the 
crossover region ranges over <j> between <f> c — (uje e r /ac) 
and <j) c + (oj£ e r /ac)- This means that, as a function of 
frequency, a sample with <fi close to 4> c will enter an ex- 
perimentally observable crossover region for sufficiently 
large u> and eventually obey Eq. (5) with a frequency 
dependence of o//( s+t ) for both the real and imaginary 
parts. 

In a series of papers reporting on the dc || and ac 
conductivities of three Graphite-Boron Nitride sys- 
tem systems it was shown theoretically that 

1. the scaling relations derived from the first order 
terms of Eq. ([!]), above and below <fi c , had the fre- 
quency dependencies required by the classical scal- 
ing ansatz @- |] (dc scaling is dealt with in §]); 
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2. the second order scaling relations, derived from 
Eq.(|l]), in the crossover region, had the frequency 
and ((j) — (p c ) dependencies required by the classical 
scaling ansatz g|- [|J; 

3. the second order scaling relations, in the region 
where the first order terms obey Eqs.(2) and (3), 
were not in agreement with the classical scaling 
ansatz [[|; 

4. the loss or conductivity term in the dielectric state 
had a frequency exponent of (1 + 1) jt and not 2 as 
required by the classical ansatz j2|- Q. 

Experimentally it was shown that 

1. the first order experimental results for the G-BN 
systems could be scaled onto separate curves above 
and below 4> c and that these scaled curves could be 
fitted by the scaling functions generated by Eq. ([!]) , 
for all (j) and u> values || . This was also found to be 
the case by Chiteme and McLachlan Q| for "cellu- 
lar" percolation systems, some of the dc properties 
of which are given in JTo| ; 

2. after taking the loss (conductivity) term of the di- 
electric component into account the exponent for 
the loss was much closer to (1 + t)/t than to 2 ||; 

3. the real dielectric constant, when plotted as a func- 
tion of <f> for the fixed lu, does not peak at 4>c but 
continues to increase reaching a peak at a value 
of <f> greater than <f> c . The resultant curve has a 
hump like appearance and the position of the peak 
depends on u>, and ac ■ It was then shown |f| 
that this behavior could be qualitatively modeled 
by Eq.@. 

In the present paper progress is reported on this last 
point. We begin with a further theoretical investigation 
of the behavior of the dielectric constant e just above 4> c . 
Since <7m is in general the result of finding numerically 
the root of a transcendental equation, this can only be 
achieved by expanding a^i around <p c . With the ansatz 



l/t 



'mc 



6 we obtain from Eq.(fiT) 



i/t 1/ 

T MC a I 



'(A 



1) + A^4 / V^ 



A{t/t 



l/s 



u<- -t/sA<t>ojr rojjfc 1/1 - (A - 1)<7j 



l/s 



(6) 



with A<f> ~ (</)— tp c )/(j> c . Note that 8 does not vanish when 
A<p = as we have to take into account additional terms 
of lower order which do not vanish for <fi = <ft c . These 
terms are omitted in Eq.(^) but they are of the same 
order as terms linear in cf> — <j) c and must be incorporated 
for reasons of consistency. Inserting the right hand side of 
Eq.(Q) for omc it is now straightforward, albeit tedious, 
to obtain the position of the maximum of the imaginary 
part. The calculation yields for the maximum 
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We emphasize that the deviation of max from <p c 
is obtained only when the expansion in powers of 
(c r //(7c) 1 / ( - s+t ' 1 is taken beyond the first power; in fact, 
the result given can be obtained only if the expressions 
are expanded up to the third power. 
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FIG. 1. Comparison of the analytic expression gm {<t>) with 
the numerical solution for 10 Hz (left) and 1 Mhz (right). The 
curve for the analytic ctm(0) has been drawn only for values 
of <f) where the analytic expression for 8(<f>) is supposed to be 
reliable. The parameters used are given in the text. 

We have checked the reliability of this analytic ex- 
pression by comparing with numerical solutions over a 
wide range of frequencies and are satisfied with its perfor- 
mance (see Figs.l). The parameters chosen in the illus- 
tration adhere to the Niobium Carbide system || which 
exhibits pronounced maxima beyond <j) c = 0.065; they 
are t = 5.46(4.71), s = 0.75(0.40), a c = 1.07- IO^to)- 1 
and <7j = — iLue e r with e, r = 5.4(6.7) for the left (right) 
case. The difference between the two curves, even at 
4> = 0c, reflects the error made by considering linear 
terms in 6 only. Nevertheless, the analytic expression al- 
lows important conclusions to be drawn for the position 
of the maximum. Note that ^> ma x does not depend on in- 
dividual values of <tj or ac , in other words the maximum 
position depends only on the ratio 07 /ac (or equivalently 
on LJEoe r /ac)- Note in particular that 

1. the deviation of ^ ma x from <p c starts with the second 
order term in (07 /crc) 1 ^ s+t ' ) ; recall that this first 
order term determines the width of the cross over 
region @- §; 

2. the larger the ratio aj/ac the further is the maxi- 
mum pushed away from the transition point (j) c . In 
turn, the position of the maximum tends towards 
4> c for a 1 j ac — > as it should; 

3. for finite values of the ratio aj/ac the distance 
^max — 4>c increases the more rapidly the more pro- 
nounced the inequality s + t > 2; 

4. to lowest order, the frequency dependence of of the 
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distance </> ma x — <fic is given by u> s +* . 
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The derivative of SJctm at <p c is easier to obtain, but the 
result is slightly more involved. We here report the es- 
sential result 
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We stress that the derivatives of SRom and Ssctm at (f) c 
are always continuous and tend to zero for {aj /ac) ~ * 0. 
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FIG. 2. Fits of experimental data at 1 Hz, 10 Hz, 1 KHz 
and 1 MHz (top left to bottom right) for conducting FezOn 
grains in a wax coated talcum powder matrix. The parame- 
ters used are given in the text. 

As it has now been established and confirmed exper- 
imentally in the present paper and in ||, Eq.(|l]) gives 
rise to the hump in the dielectric constant just above 4> c . 
However, as can be seen in ||, Eq. ([[]), as it stands, only 
qualitatively models the experimental data. In order to 
better fit the dielectric data we propose an effective de- 
pendence of ac on (j) that is based in spirit on a model 
due to Balberg |ll[ for non-universal values of t. The 
model accounts for values of t higher than those allowed 
by the Random Void (RV) and Inverse Random Void 
model (IRV) |L2|,[L3). In this model Balberg assumes that 
the resistance distribution function /i(e), where e is the 
proximity parameter, has the form e~ w as e — > and 
does not tend towards a constant as in the RV and IRV 
models. Using the approach of Jl3],[l4| and keeping the 
underlying node links and blobs model |^), Balberg de- 
rives an expression for a non-universal t which is equal 
to the one obtained from the RV model when w = 
but can give a larger t for w > 0. For w > the model 
also shows that the average resistance in the network can 
diverge when <fi — > <f> c . The increase in t beyond its uni- 
versal value t un as found by Balberg can be lumped into 
t = t un + t nun + r where r is the extra contribution due 
to the characteristic resistance of the network diverging 
at <fi — > <f) c . This increase in characteristic resistance is 
incorporated into Eq.([l]) and hence also into Eqs.(||) and 
(||), by substituting ac with 



C — a 00 + °c0 



with r > 0. The solution of Eq.(|l|) should yield a contin- 
uous a \i across <p c , which means that the conductivity 
a c should nowhere vanish. As a consequence, ac must 
reach a non-zero value denoted by (Too at (j) c which is 
expected to be considerably lower than a C Q. To avoid 
too many parameters the simplest assumption was made 
for (ft < <p c , that is ac = coo This should be reasonably 
valid just below <fi c where our interest is focussed. We 
stress that the effect of the modification by Eq. (||) is vir- 
tually indiscernible with regard to first order effects for 
the solution om of Eq. (nh or the corresponding percola- 
tion power laws (Eqs. (2) and (3)) in the region where 
the power laws are obeyed. Computer simulations show 
that using Eq.(||) or a constant ac in Eq.(|l|) causes a 
difference in a M (4>-4>c) {4> > 4>c) or e M {4>c - 4>) (4>c > 4>) 
which is too minute to be resolved from available experi- 
mental data. The same holds for dispersion plots against 
lu of cm above (f> c and £m below 4> c as given in [Q. For 
this reason it has never been necessary previously to con- 
sider the modification given by Eq.(||). However, in our 
context the higher order effects are discernible in the di- 
electric constant just beyond <f> c and can be fitted much 
more satisfactorily than in Q using Eqs.((l|) and (|). 
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FIG. 3. Fits of experimental data at 1 Hz for T = 25° 
(top left) and T = 120° (top right). The bottom row displays 
corresponding results for 1 KHz. The system is the same as 
in Fig. 2. The parameter values are given in the text. 

In Fig. 2 we display data and corresponding fits for con- 
ducting Fe^O^ grains in a wax coated talcum powder ma- 
trix J0jn]] ■ The fits are far superior than those presented 
in pjT We have chosen (somewhat arbitrarily) (Too = 
crc/lO with the values ac = a c o = 2.63 ■ lCP^fim) -1 
obtained from extrapolating dc experimental results to 
0=1, and a j = —iu)e§e T where e r is measured sep- 
arately at each frequency. Good fits are obtained for 
(t-r,r,s) = (4.7, 0.5, 0.97), (4.3, 0.7, 0.98), (4.4, 0.6, 1.0) 
and (5.6,0.4,1.6) when moving from the top left to the 
bottom right display in Fig. 2. The value 4> c — 0.025, 
obtained from dc measurements, is used in Figs. 2 and 3. 



Note that t = t UD + t n 



r is somewhat larger than the 
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separately measured non-universal dc value of 4.2. Also 
note that in all cases it turned out that r < 1. This 
implies a steep increase of a^ just above <f> c . The fits 
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are not necessarily optimal as the multi-parameter land- 
scape of the least square expression |cr^ p — <r^ m | 2 in the 
parameters t, s and r (for fixed (Too) has many local min- 
ima. However, the quality of the fits at the different local 
minima does not vary greatly for acceptable fits. 

While this second order behavior of e is in principle a 
complicated function of the various parameters the po- 
sition of the maximum at max is essentially dependent 
only on the quotient ai/ac- This implies that a decrease 
of frequency (which is a decrease of 07) has an effect 
similar to a corresponding increase of ac, which can be 
obtained by an increase of temperature. This is convinc- 
ingly demonstrated in Fig. 3 where fits are obtained for 
the same system as in Fig. 2 at different temperatures, i.e. 
at different values of ac and e r . The experimental data 
and the theoretical curves show that an increase in ac 
is equivalent to a decrease in to (or 07). The experimen- 
tal data appear somewhat erratic as they usually do for 
the low frequency chosen (1 Hz). The data in the first 
column of Fig. 3 arc slightly different from those in the 
first column of Fig. 2 as the former are taken in a tem- 
perature controlled oven. Much better fits are obtained 
for higher frequencies as illustrated in the bottom row of 
Fig. 3 (1 KHz) as the dielectric data are more reliable at 
higher frequencies. Similar to Fig. 2, experimental values 
are used for e r (u>, T) with a c o — ac = 3.22 ■ 10 _1 (f2m) _1 
at 25°C and a c0 = a c = f.42(^m)- 1 at f20°C. The 
fits were obtained with the parameters (t — r, r, s) — 
(4.65, 0.5, 1.0), (3.4, 0.6, 1.0), (4.4, 0.6, f .0), (3.5, 0.9, 1.0) 
again choosing the fixed value (Too = fc/10. 
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FIG. 4. Fits of experimental data at 100 KHz of a wa- 
ter-oil emulsion using a hydrophil agent for uniform droplet 
size. 

In continuum systems, all models for a non univer- 
sal t are based on the premise 0- @ that a distri- 
bution of conductances exists in the conducting compo- 
nent which has a power law singularity. In a system 
where no such singularity occurs it should be possible 
to fit the data using Eq.(l), with t values closer to 2 
and without the use of Eq.(8). A system par excellence 
which satisfies this requirement is a water-oil emulsion, 
using a surfactant to ensure uniform drop size. Data 
for such systems has been published in [16] -[18] and we 
present in Fig. 4 data from Fig.l of [16] and from Fig. 4 of 
[17], which have been fitted using Eq.(l) with the same 
value of ac above and below <j) c . The values used, i.e. 
a c = 1-9(1. 5)(fim)-\ 07 = -i2 ■ HT 11 ^ (f> c = 0.2(0.06) 
for the left (right) display and v = 10 5 Hz are based 
on the dielectric and conductivity curves given in the 



quoted papers. The fitting curves use s — 1.15(1.35) and 
t = 1.85(2.1). Note that the values of t are close to the 
universal value. Although better fits should be obtained 
with a more precise knowledge of the system parameters, 
these results show that using Eq.(8) as input in Eq.(l) is 
necessary to fit the dielectric data only when abnormally 
high values of t are observed in the dc data. Note that 
these data have been previously analyzed using Eq.(l) 
with a single exponent (s = t) in [19]. 

To summarize: the pronounced hump observed exper- 
imentally for the dielectric constant of percolative metal 
insulator systems just above the critical concentration 
can be modeled using the higher order terms of Eq.(l) 
near <j> c . In the cases where t is well above the universal 
value Eq.(8) must be used to quantitatively model the 
dielectric data. Note also that all previous experiments 
on the dielectric constant "below" <j) c , where <fi c has not 
been independently measured by dc conductivity mea- 
surements, have probably incorrectly identified values of 

0c- 
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